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1.0 INTRODUCTION 

Incorporation of the four gimbal I MU model into the IMUCAL bench program 
has Isolated an anomaly In the star tracker software line of sight (LOS) 
rate test. During verification testing of the revised IMU model, It was 
noted that the LOS rate cosine was greater than one In certain cases. 

This implies that either one or both of the star tracker measured end- 
point unit vectors that are used to compute the LOS rate cosine have lengths 
greater than unity. The search for the software rogue that is stretching 
these vectors came to an end at the roll/pitch nonorthogonality matrix 
in the TNB CL module of the IMU software. 
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2.0 DISCUSSION AND RESULTS 

The star tracker software (Ref. 1) generates three measurement vectors 
for each star that Is tracked. Each vector Is rotated through the navigation 
base to inertial platform transformation matrix that is computed In the 
TNB_CL module. One of the three vectors is computed from averaged rueasurement 
data and Is used for IMU alignments. The retraining two vector?, are computed 
from the first and last measurement samples jf the twenty-one sample set 
and are used Internally in the LOS rate test. Any corruption of the vector 
to be used for alignment is of primary importance. 

The source of the anomaly is the roll/pitch gimbal nonorthogonality (DP) 
matrix in the TNBjCL module (Ref. 2). The ideal matrix 


cos5 ~$in$ 0 

(E) * sinfi cos-$ 0 
t 0 0 1 


is approximated in the current software design by 



where 6 is the roll/pitch gimbal nonorth, jonality, DP. In other words, 
the nonorthogonality matrix is n.northogonal. 


2.1 EFFECT OF THE DP MATRIX APPROXIMATION ON THE IMU ALIGNMENT ACCURACY 

Star measurement vectors that are used for IMU alignments are unitized 
in the aberration correction equations and, therefore, are not corrupted 
in length. 

However, there is some concern that the DP matrix may introduce error 
into the measurement vector direction. If e is defined to be the angular 
error in the measurement vector direction, then 


cos e = UNIT{(A)$} • 

(Vj 2 + y 2 2 H cos5 + 5 sin5 ) + v 3 2 

cos e = — 

((1 + <$ 2 )( Vi 2 + V2 2 ) + V 3 2 )!/ 2 


The error, e, will be greatest when V 3 = 0 , therefore. 


e max * 


cos 


-1 


cos 6 + 6 sin 6 
(1 + 62 ) 1/2 _ 


(1) 


which is plotted as a function of <$ in Figure 1. Since the range of 
for the Shuttle IMUs is |$| < 0.002 radians. Figure 1 shows that measurement 
vector directional errors induced by rotation through the OP matrix are 
insignificant. 


2.2 EFFECT OF THE OP MATRIX APPROXIMATION ON THE LOS RATE TEST 

The end-point measurement vectors are not corrected for aberration and, 
therefore, are not unitized after rotation through the DP matrix. The 
LOS rate test assumes that these vectors are of unit length and computes 
the LOS rate cosine by 


cos p = il • V 


where P is the star LOS "rate", U the first measurement vector, and 
V the twenty-first measurement vector. Since the end-point vectors are 
< greater than unity in length. ^th^ computed LOS rate is greater than the 

ideal value by the factor 1/1 U | 1 V | . Measurement data are accepted whenever 


U • V > T0L12 


Therefore, the probability of accepting deDris data as star data is increased 
due to this anomaly. The end result of this error is to make the LOS 
rate test less sensitive. 

The angular error, e, corresponding to the cosine error can be computed 
by comparing the flight software formulation with the ideal solution. If 
P is defined as the corrupted inner product using the approximation of 
the DP matrix and l 5 the ideal inner product, then 
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P = (A# • (A)V 
P « (A# * (A)? 

which in expanded form are given by 

P = U1V1 + U2V2 + U3V3 + 6 2 (UiVi + U2V2) 

P = U1V1 + U2V2 + U3V3 

The error in P will be greatest when 113 = 1/3 = 0 . Assuming these worst 
case conditions (U3 = V3 = 0 ), F and P take the form 

P = 0 • $(1 + <5 2 ) = cos p 

"P = $ • $ = COS P 

Therefore 

cos p = tf • if(l + 6 2 ) = (1 + <5 2 ) cos p (2) 

The error in the equivalent angle (e * |p - o|) will be greatest when 
p a 0; however, if we set p = 0 we cannot solve for p because 

p = cos'^(l + 6 2 ) 


is undefined for 6^0. 

Alternatively, we can set the computed angle, p, equal to zero and solve 
for the actual angle, p, that would yield a computed angle of zero. 



(3) 


The maximum equivalent angular error in the computed star LOS rate is 
plotted as a function of the roll/pitch gimbal nonorthogonality angle 
in Figure 2. Also plotted are average errors, em e an» that were generated 
by a simulation program. Figure 2 illustrates that the DP matrix approximation 
introduces large errors into the computed star LOS rate. 
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The significance of the LOS rate error Is more easily understood by expressing 
the error in terms of a LOS rate sensitivity, x* IT the DP matrix approximation 
were replaced with the exact matrix, the LOS rate sensitivity would be 
defined as 



where TOL12 is the LOS rats test tolerance limit and 3.2 sec Is the time 
between the first and last star sightings. By substituting T0L12 for cos p 
in equation (22, the largest angular change between the first and last 
measurements, p, is given by 



The LOS rate sensitivity, therefore, is 



The LOS rate sensitivity is plotted as a function of 6 for three cases 

in Figure 3 (solid lines). The uppermost curve corresponds to the current 

STS-1 software configuration (the STS-1 value for T0L12 is padded to circumvent 

the unforeseen). The dashed line with the same X-axis intercept illustrates 

the LOS rate sensitivity that could be obtained by using the exact DP 

matrix or by unitizing the end-point measurement vectors before forming 

the inner product. The middle curve was generated using the minimum recommended 

value for T0L12 (Ref. 3). This case corresponds to the limit of the current 

formulation of the LOS rate test. As with the previous case, the dashed 

line indicates the capabilities that could be obtained by using the exact 

DP matrix or unitizing the end-point measurement vectors. Since the present 

inner product algorithm is inherently inaccurate, the bottom curve is 

plotted to illustrate capabilities that could be realized by using the 

vector cross product of the end-point measurement vectors to compute the 

sine of the angle change. The cross product algorithm has proved to be 

very accurate; and, in addition, its accuracy is not degraded by the approximation 

of the DP matrix. The bottom curve, therefore, defines the LOS rate sensitivity 

limit imposed by the hardware inaccuracies. 
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Returning to tlw STS-1 cast, the Ideal LOS rate sensitivity (ngn* 4) Is 


Mdeali - 0.041 deg/sec 


Tlw eaxlnwi arror resulting from tlw OR matrix approximation (5 - 0.002) 
limits tlw LOS rate sensitivity to 


X - 0.065 deg/sec 


This limit Is 596 greater than the Ideal limit and, therefore, the 10S 
rate test will accept significantly more debris as alignment data as a 
result of the errors in the OR matrix approximation. 





3.0 CONCLUSIONS AND RECOMMENDATIONS 

The approximation of the 0? Matrix does not corrupt the directional Information 
In star measurement vectors and* therefore, will not degrade IHU alignment 
ccuracy. The OP matrix approximation does, however. Increase the length 
of vectors that are rotated through It. As a result, significant errors 
are Introduced Into the star LOS rate test. Therefore, the possibility 
of accepting debris position vectors as measurement data for I MU alignments 
Is Increased. It Is recomended that the approximation of the OP matrix 
in the TNBJCL module be replaced with the exact matrix 

f CDP - SOP Ol 

SDP COP 0 


where SDP and CDP are the sine and cosine, respectively of the outer roll/ 

E itch gimbal nonorthogonality angle. SOP and COP are currently supplied 
y the IMU ground calibration software and, therefore, additional flight 
software code Is not required to compute these constants. The nine equations 
for the elements of the Nav base to cluster transformation matrix, however, 
would have to be rederived to include the gimbal nonorthogonality sines 
and cosines, and then these equations implemented into the flight software. 

In addition to this reconnendatlon, two additional solutions exist. Unitizing 
the end-point measurement vectors would restore their lengths to unity 
thereby eliminating that source of error. The second alternative is to 
reformulate the star LOS rate test using the cross product algorithm. 

Choice of one of these alternatives will be made via discussions with 
RI, JSC and IBM in the near future. 
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